EVOLVING COMMUNITIES WITH INDIVIDUAL PREFERENCES 
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Abstract. The goal of this paper is to provide mathematically rigorous tools for modelling the 
evolution of a community of interacting individuals. We model the population by a measure 
space (QjJ 7 , v) where v determines the abundance of individual preferences. The preferences of 
an individual u> £ Q are described by a measurable choice X (uj) of a rough path. 

We aim to identify, for each individual, a choice for the forward evolution Yt (ui) for an individ- 
ual in the community. These choices Yt (uj) must be consistent so that Yt (ui) correctly accounts 
for the individual's preference and correctly models their interaction with the aggregate behaviour 
of the community. 

We focus on the case of weakly interacting systems, where we are able to exhibit the existence 
and uniqueness of consistent solutions. 

In general, solutions arc continuum of interacting threads analogous to the huge number of 
individual atomic trajectories that together make up the motion of a fluid. The evolution of the 
population need not be governed by any over-arching PDE. Although one can match the standard 
nonlinear parabolic PDEs of McKean-Vlasov type with specific examples of communities in this 
case. The bulk behaviour of the evolving population provides a solution to the PDE. 

An important technical result is continuity of the behaviour of the system with respect to 
changes in the measure v assigning weight to individuals. Replacing the deterministic v with the 
empirical distribution of an i.i.d. sample from v leads to many standard models, and applying 
the continuity result allows easy proofs for propagation of chaos. 

The rigorous underpinning presented here leads to uncomplicated models which have wide ap- 
plicability in both the physical and social sciences. We make no presumption that the macroscopic 
dynamics are modelled by a PDE. 

This work builds on the fine probability literature considering the limit behaviour for systems 
where a large no of particles are interacting with independent preferences; there is also work on 
continuum models with preferences described by a semi-martingale measure. We mention some 
of the key papers. 

Key words and phrases : Rough paths analysis, interacting particle systems, propagation of 
chaos, dynamic economic equilibrium. 

1. Introduction 

Consider a community = {ajj , i e 1, . . . , N} of N individuals who at time zero are allocated 
positions (Yq (w)) IiIg Q m some state space, and suppose that these individuals have preferences which 
determine how they evolve in their environment. Let the evolution of individual Ui be denoted Y 1 , 
and suppose X % is the preference of this individual. To model this situation we assume that Y l and 
X 1 are related via 

dY l = f (Y*) dX\ 
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This equation has a unique meaning when X 1 is a p— rough path and / is Lip-7, 7 > p (see |36j). 
By imposing a measure on ft (e.g. the counting measure), we can introduce T =Law(Y), and then 
consider the situation where the evolution Y l is influenced by the wider population through T. This 
leads to equations of the form 

dY = f (Y) dX + g (Y, V) dt. 

Deterministic models of this type are commonplace in the modelling of physical systems. For 
example, our individuals might be planets. Although planets do not have preferences, they are 
subject to gravitational forces exerted by other planets which affect the state Y = (p, q) of their 
position and momentum. In this setting, we can model the Newtonian evolution of the locations of 
the individuals via a differential equation: 

apt (uj) = - > Gm (w) m (w) - — — —3 dt 

zen\{u} \Qt(u)-qt(u)\ 

dq t (uj) = — r^Pt ( w ) dt 
in (u>) 

G = 6.6730010~ n N(m/kg) 
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Unless there is a collision, the theory of ordinary differential equations is an adequate tool to 
describe the short time evolution of this system. If each member of the community carried charge 
of the same sign then the equations would change again 

dp t (uj) = V . I — C{u)C{u) qt{uj) ~ qt{[d \ p t (Lo)xd qt (£o) 

dq t (u) = — T^Pt ( w ) dt 
m [ui) 

C is the charge on to and /io=47rlO _7 N/A is the magnetic permeability. 

These examples are very different in detail because of the involvement of the dqt ■ But both equations 
capture systems of physical interest, where it is natural to consider the evolution of the population 
as a whole, and to understand what happens when particles are replaced by more particles with 
proportionately smaller mass (or charge) in the continuum limit. In this case one would hope and 
expect that the particle and/or current density would solve the appropriate Vlasov style equation. 

Of course, there are a huge number of similar if less precisely characterised models: in the social 
sciences, in the modelling of the evolution of cancer, etc. where the evolution of an individual is 
affected by the evolution of the wider community. It may not always be the case that the interaction 
with the population is pairwise and more generally, one might expect to consider equations of the 
form 

dY t (u) - (Y t (w) , v*Y t ) dt, 

where v*Y t is the push forward of v giving the mass distribution of Y at time t, and Y t (uj) is a 
solution to the above equation for /i-every uj. 

Individual differences mean that different individuals will respond differently to the same external 
environment. We can easily make adaptations to the calculus to take such behaviour into account. 
One is lead to equations of the following kind: 

dY t (oj) = <j> {Y t (w) , v.Y t )dt + £ (Y t (w) , u*Y t ) dX t (w) 



EVOLVING COMMUNITIES 



3 



where X (oj) represents the individual preferences of the individual uj. Now if X is smooth there is 
no additional difficulty. If v is a probability measure, and X is a semi-martingale under this measure 
then (under regularity conditions) Sznitman [42] , Kurtz [27] , Meleard [37] and others have proved 
that the corresponding particle system obtained by taking an i.i.d. sample from X and using this 
empirical measure in the above equations provides a converging sequence of particle systems. The 
limit can be identified with the law of a non-linear PDE which solves the Vlasov equation. Dawson 
and Gartner [9] have important results on the large deviations in the convergence of the weakly 
interactive system (where £ (Y t (ui) , v^Y t )=^ (Y t (w)) and later den Hollander [7] and Guionnet [T], 
|llj.|12j and |23) considered the large deviations for interaction in a random media in problems 
arising from the dynamics of spin glasses. Kurtz promoted more advanced discussion in [28] . 

In many cases of interest, it is unreasonable to expect the preferences X to be a semi-martingale 
as evidenced by the sucess of fractional Brownian motion in the modelling of fluids (see [24] and 
the references therein). In addition, individuals often have knowledge that makes the previsible 
assumption equally inappropriate. We now understand that the natural assumption on X that 
leads to equation with a strong meaning is that X should be rough path. Indeed there are a 
large number of deterministic (and numerically approximablc) systems that evolve without the 
assistance of a PDE. We will study a mathematical framework which exposes the consequences of 
persistent differences between individuals in the population dynamics (see [41] for a study of such 
a phenomenon in the context of red deer populations.) 

The McKean- Vlasov model leads one, in the limit, to the equation 

dY t (w) = (Y t (w) , u.Y t ) dt + Z (Y t («) , u.Y t ) dW t (u) 

where the individual preferences are given by a d-dimensional Wiener measure W. However indi- 
viduals can have very different volatility and speed of reaction to events. Let a be a positive real 
function on the space fi of individuals and consider the equation 

dY t (u) = <t> (Y t (uj) , i/,r t ) dt + i (Y t (lo) , i/,y t ) a (w) dW t (uj) . 

For appropriate <j), £ and paths fi s in measures on K-space one can consider the indexed family of 
differential equations, one for each lo, 

dY t (w) = (Y t (u) ,iH)dt + £ (Y t (w) , Ht) o M dW t (lo) . 
Y (lo) given 

For almost every ui the path t —¥ a (lo) Wt (ui) is a geometric rough path of finite p-variation for 
every p > 2 based on rescaling W and its Levy area. If v t is a path of finite variation in the space 
of measures and <f>, £ are at least C 2+e then it will be the case that the rough path solution Y t (uj) 
to this equation will exist and be unique. Considering all a;, we see that Y t (uj) is a random variable 
and we denote its law by the probability measure p t . Of course, this new path t —> p t in measures 
will not in general coincide with the path t — > fit ■ But it makes complete sense to ask whether there 
is a choice t — » pt so that the resultant measure path t — > jit does coincide with it. In this case we 
have a community of individuals evolving according to their individual preferences in a way that is 
also consistent with the dynamics of the population as a whole. 

We note that in general having individuals with different volatility results in a process t — > 
a (uj) Wt (u>) that is far from a semimartingale against the Wiener measure and using the base 
filtration. One cannot have a (uj) measurable in J-q unless one enlarges the filtration or a is constant. 
The lack of prcvisibility docs not impede the rough path perspective, and there is no issue at all in 
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setting up the equations. One theoretically amusing choice for a is to take 



su Pte[o,i] \W t (u)V 

which in some sense eliminates enthusiastic outliers in the population. 

To move on from posing a meaningful questions to identifying solutions is actually quite chal- 
lenging. For example it is not clear, at the level of generality that we introduce, that the path 
t jit will have bounded variation or what Banach space to consider it as a path in even if it does. 
This raises another issue - in that solving equations such as this we require the pair (/x*, Wt (ui)) to 
be a rough path which normally requires extra data unless one has good control on t —> fj,t so we 
have to make some compromises. No doubt there is much that can be refined and taken further. 

Let t —¥ [it be a path in the space of probability measures representing a putative evolution of 
the population Y t (oS) . We introduce the "occupation" measure process T t := L /i t dt and note that 
it is monotone increasing and Lipschitz with norm one in the total variation norm on measures. 
Let 9 (y) [i — J 4> (y, y') M {dy') , then 9 can be viewed as a linear map from our space of measures 
to vector fields on the Y— space. 

We make two significant simplifications to make the problem more tractable: 

(1) We only allow so-called weak interactions between the individual and the population which 
take place only in the drift component of the equation i.e. £ (Y t (u>) , v*Yt)=£ (Y t (u>)) . 

(2) The interaction between the individual and the population admits a superposition principle. 

Together these imply that we can write the interaction between Y, its preferences W and the 
distribution v of the community in the following form 

dY t (w) = 9 (Y t (w)) dT t + £ (Y t (w)) dX t (w) . 

We then look for fixed points of the map that takes t — > /it to t — > jit . Since T has bounded variation 
this equation poses fewer technical problems than the general case but still allows discussion of the 
Vlasov type problems discussed initially. 

The paper is structures as follows. In Section[2]we spend some time setting up our notation for the 
rough path framework; this is the mathematical technology we use to model the community. Section 
[3] then explores the special case where the law of the preferences is given by a finitely-supported 
probability measure on the space of (geometric) rough paths. Here we prove an existence and 
uniqueness theorem for the law of the nonlinear McKcan- Vlasov RDE. The methodology here is 
distinct from that used later to prove the (more general) result for non-discrete measures. But this 
simple case allows us to see very clearly how the weak interaction assumption, combined with the 
LV Extension Theorem of [33] gives rise to the uniqueness of fixed points. In Section [4] we proceed 
with the roadmap sketched out above. We first present some Gronwall inequalities for rough 
differential equations, developing the deterministic estimates from [20j and focusing particularly 
on the conditions needed to ensure the integrability of the estimates. We make use of the recent 
paper [5] in showing that these conditions arc satisfied for a wide range of preference measures. We 
then present conditions that ensure the existence and uniqueness of fixed points, and discuss the 
their continuity in the measure on preferences. Finally, Section [5] establishes propagation of chaos 
(a la Sznitman [42]) for the convergence of the finite particle system. We note that this paper has 
already lead to follow-up work (see, e.g., [2]); we discuss other possible applications of our results. 
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2. Preliminaries on rough path theory 

There are now a wealth of resources on rough path theory, e.g. [32], [20] , [HJ , [36] . Rather than 
give an overview, we will focus on the notation we need for the current application and direct the 
reader to references where appropriate. We first recall the notion of the truncated signature of a 
parameterised path in C 1 ~ var ([0,T] , R d ) (the set of continuous paths of bounded variation), this 
is: 

N 

S N {x) s t := 1 + ^ / dx ti ® dx t2 <S> ... ® dx tk € T N (R d ) . 

fc=l J 9<tl<t 2 <. — <t k <t 

Where T N (R d ) = © ^ (R d ) ^'denotes the truncated tensor algebra. We use 7r„ to denote the 
canonical projection 

7T n : T N (R d ) -> n = 0,l,...,iV. 

For x™ in (R d )® n we define x n '(* lv "'* n ) to be the real number 

x »,-(u,..A0 = ... g, e *J (x nj =; (e^.^.x") , 

where e*, denote the standard dual basis vectors. We equip each (R d )®" with a compatible 
tensor norm |-|( K d)«™, and let 

c?Ar(g,h):= max \m (g - h)| fRd -,®i . 

2= 1,...,7V v ' 

It is a well-known that the path Sn {%) in fact takes values in the step-TV free nilpotent group 
with d generators, which we denote G N (R d ) . Motivated by this, we may consider the set of such 
group- valued paths 

x t =(l,x t V..,xW)eGW (R d ), 

for p > l.We can then describe the set of "norms" on GLpJ (R d ) which are homogeneous with 
respect to the natural scaling operation on the tensor algebra (see [20] for definitions and details). 
The subset of these so-called homogeneous norms which are symmetric and sub- additive (|20j) 
give rise to genuine metrics on GLpJ (R d ) . And these metrics in turn give rise to the notion of 
a homogeneous p-variation metric d p . vax on the GL p J (R d )-valued paths, a typical example being 
the Carnot-Caratheodory (CC) metric dec- The group structure provides a natural notion of 
increment, namely x Sjt := xj 1 ® x t and we may then define 

Vvar;[o,T](x,y):=||x-y|[ [0T] := sup <fcc(xt i ,t, +1 ,yt,,t, +1 ) ■ (2.1) 

If (j2~Tj) is finite then, w x (s,t) := | |x| \p_ var . [s t] (:= ||x- l||£_ var .[ S)t ]) is a controfl. Also of interest 
will be the inhomogeneous rough path metric defined by 

(\ i/p 

M.e. it is a continuous, non-negative, super-additive function on the simplex A[q T j = (s, t) : = s <= t = T 
which vanishes on the diagonal. 
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And the uj— modulus inhomogeneous metric, with respect to a fixed control oj, which is defined by 

/ \ i i \n (x a ,t -y s ,t)\ 

Pp-u,-,[o,T\ (x,y) = |x - yo| T Lpj (R< i) + max — . 

t=i,...,Lpj uj[s,t) ,p 

The space of weakly geometric p— rough paths will be denoted WGVL p (R d ) . This is the set of 
continuous paths with values in GLpJ (parametrised over some, usually implicit, time interval) 
such that (|2.ip is finite. A refinement of this notion is the space of geometric p— rough paths, 
denoted Gft p (R d ) , which is the closure of 



{S M (x) o y-x&C 1 - mr ([0 7 T},R d )} 



with respect to the rough path metric d p _ var . 

We will often end up considering an RDE driven by a path x in WGil p (M d ) along a collection 
of vector fields V — (V 1 , V d ) on R e . And from the point of view of existence and uniqueness 
results, the appropriate way to measure the regularity of the ViS results turns out to be the notion 
of Lipschitz-7 (or, simply, Lip-7) in the sense of Stein. This notion provides a norm on the space 
of such vector fields, which we denote \-\ Lip _ 1 ■ We will often make use of the shorthand 

\V\ Lip ^=.m^ d \Vi\ Lip _ 7 . 

Finally, throughout the article we will consider spaces of probabilities measure on metric spaces 
(S,d). 

Notation 1. We will use A4 (S) to denote the space of probability measures on (S,B (S)) . For 
p > 0, AA P (S) will represent the subset of hA (S) which have finite p th -moment in the sense that 



d (so, s) p /1 (ds) < 00, 

JS 

for some (and hence every) Sq G S 

It will be convenient to have a shorthand notation for some of these spaces. 
Notation 2. We will write 

(S Nte ,a N ) for (G N (R e ) , d N ) , and (P p , e , p p ) for (Gfl p (M. d ) , p p - va r;[o,T]) 
Furthermore (St, e ,&t) will mean (iSui e , crujj whenever t is not an integer. 

3. Weakly interacting communties 



Let (Mt) te r t] b e a family of probability measures on GLpJ parameterised by time. The main 
object of study in this paper will be solutions to rough differential equations which incorporate weak 
mean-field interactions with (fit) te [o t] ■ By this we mean equations of the following type 

dy t = a(y t ,iT 1 y)iJ l t(dy)dt + V(yt)dx,y(Q) = yo. (3.1) 

The rough path x flows along the vector fields V = (V 1 , V d ) , but the resulting trajectory is also 
influenced by (/xt) tg [ x] through the interaction kernel a. Assuming enough regularity on the path 
!h}(i we may define the integral 
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a continuous bounded variation path in an appropriately chosen ambient Banach space. It is 
convenient to rewrite the main equation (|3.1[) as 

dy t = V° ( Vt ) dr^ + V (Vt) dx t + V (y t ) dx t , y (0) = y , (3.2) 

where V° and a are related by 

V°(y)(fi) = f *(y,ir iy )[x(dy). 

We will discuss the detail of this construction in Section[U In the cases we consider, (/J.t)te[o t] wm " 
be derived from the marginal distributions of a probability measure in M. (P Pt e) ; i-e. those derived 
from pushing- forward under the evaluation maps ipt (x) = x t , t £ [0,T] . We denote a solution to 
(|3.2[) by Qv°,v (Mj Uo, x ) , and fix a probability measure uj x j/ on K c x GT2 p (lR d ) . uo describes the 
initial configuration of the particles and v is the law of the preferences or, more conveniently, the 
preference measure. By taking a realisation (Yo,X) of uo x v on some probability space (f2, J 7 , P) , 
and then using X to solve (|3.2j) we will have constructed a well-defined map ^ from the space 
M (Gflp (K e )) to itself given by the push-forward: 

: [L h-» [9yo,y (/i, •, •)] „ ("o x f ) • 

will then be fixed point of this map if and only if O^oy (/x, IojX) is a solution the (nonlinear) 
McKean-Vlasov-type RDE 

dY t = V(Yf)dX t + V° 

Law (Y) = /i, Law(y ) = u o ' 

A key objective of this paper is to demonstrate that there exist unique fixed points to (|3.3[) for a 
class of preference measures which extend far beyond the usual semimartingale setting. 

We first spend time developing an important special case, namely when v is a finitely-supported 
discrete measure of the form 

N 
i=l 

In this setting, we can attempt to resolve the fixcd-point-problcm (|3.3[) by solving the system of 
RDEs 

N 

<fyi{t) = V(y i {t))(b^(t) + ^^a(vi{t),v j [t))dryjt, y<(0)= tt (3.4) 
for i = 1, JV. And then defining the measure to be the convolution 

where uf n is the n-fold product measure of uq. More precisely this means that 

N 



where we have written y^* to emphasise the dependence of on its starting point With /i 
defined in this way we would expect that ^ v (fi) = [i, and indeed this approach will work for 
smooth preferences. In the rough case (p > 2) however things are more complex. Here in order to 
solve (|3.4[) we need to define a priori the cross-iterated integrals between the (components of) the 
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preferences x 1 and x J . The LV Extension Theorem ([33]) guarantees that this can always be done, 
but in general there are many choices for the extension. To ensure uniqueness of the fixed point, 
we need to check that the resulting solution is not sensitive to this choice; the remainder of this 
section will present conditions which will guarantee this. 

The results of this section will later be subsumed by the general fixed point theorem of Section 
2J Nonetheless they are important for three reasons. Firstly they expose, in an original and lucid 
way, the importance of the weakly interacting structure; secondly, they highlight the main obstacle 
in extending the analysis to general interactions, in a way that cannot be easily discerned from 
the general fixed point result; thirdly, they crucially underlie our later treatment of the convergent 
behaviour of the finite particle system. 

3.1. A two-particle system. To make clear the structure of the argument, we first deal with the 
case where N = 2 and p G (2, 3) ; i.e. the preference measure is supported on only two geometric 
rough paths in G£l p (R d ) . We write v = XS Xl + (1 — X)S X . 2 . By the LV Extension theorem there 
exists an element x in WGVL P (R 2d ) which lifts (xi,x 2 ) consistently with xi and X2 in the sense 
that, 

Pjx 1 = xj and {Pj ® Pj) x 2 = x 2 for j = 1,2 

where P\,P<2 : K 2e = R e x R e — > M e arc defined by P\z = x and P 2 z = y when z = (x,y) , and 
where (Pj ® Pj) (zi <g> z 2 ) = PjZi <g> PjZ 2 G (K 2e ) . We can simplify this by writing 



x 1 = (xj.xj) e R. 2d , x 2 = ( *i ) G (R^f 2 , (3.5) 



under the obvious identifications. The only constraint on the terms (*) arises from the need to 
make x weakly geometric. Given such an extension, we can solve the following RDE uniquely 

dy t = W° (y t ) dt + W (y t ) dx t , y (0) = ( yi (0) , y x (0)) G M 2e . (3.6) 

Wherein W = (W 1 ,..., W 2d ) is the collection of vector fields on R 2e ~K e xR e defined by 

W* ( yi ,y 2 ) = (V* ( Vl ) , Or.)* , i = 1, d (3.7) 
W l ( Vl ,y 2 ) = (0 K e, V* (y 2 ))* , i = d + 1, 2d, 
and the interaction is transmitted through 

W" (2/1,2/2) = A (cr (2/1, 2/1) ,<r (2/2,2/1))* + (! - A ) (o- (2/1, 2/2) ,0" (2/2, 2/2))* ■ 
By writing the solution y in terms of its projections 

y 1 = (yi.yj) e M 2e , y 2 = ( y j * 2 ) e (R 2e ) 82 , (3.8) 

we can obtain y t = {l,y\,y$) G WGfl p (R e ) . We will prove that the probability measure 

It = w® 2 * [X8 yi + (1 - A) d y2 ] G M (P P , e ) (3.9) 



is a fixed point of the map We will then show that every fixed point has the form (|3.9[) ; i.e. its 
suppport is {yi,y2} , where yi,y2 arc projections of the solution to (|3.6[) driven by any extension 
x. The uniqueness of the fixed point will follow by proving that the projections yi and y2 do not 
depend on the extension (and hence neither does the measure (|3.9|0 . This is the essential content 
of the following proposition. 
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Proposition 3.1. Let 2 < p < 3 and y\ (0) ,yi (0) G R e . Suppose that Xi and X2 are iwo elements 
ofGflp (R d ) . Assume that W° and W = (W 1 , W 2d ) are, respectively, vector fields in Lip 13 (R 2e ) 
and Lip 1 (R 2e ) for some f3 > 1 and 7 > p . Let x 6e any element of WGVl p (K 2d ) which extends Xi 
and X2 in the sense of H3.5\) , and let y be the unique solution in WGfl p (R. 2e ) in to the RDE \3. 6)) 
driven by x. Then y has the property that its projections yi , y2 ( as given in \3. 8\) ) are elements of 
WGQq (R e ) which depend on xi and X2, but not on the extension x. 

Proof. We prove that yi and y 2 , the projections of the solution to (|3.6[) . depend only on Xi and 
X2 and not the iterated integral between them. In other words, that yi and y2 have meaning 
independently of the terms * needed to specify the joint lift in (|3.5[) . To see this we recall (|14j) 
that y Sj t, the increment of the path level solution over [s, t], is equal to 

mft E [ W ° - **) + W (»*«K,t j+1 + DW (y u )W (2/ tl )x 2 >ti+i ] , (3.10) 

i:ti£D[s,t] 

where DW (y ti ) W (y ti ) [x<8y] = DW (y ti ) [W (y ti ) x] [y] , and D [s, t] denotes a partition of [s, t] . 
The last term in the summands in (|3.10p equals 

\DW{y u )W{y u ) [x^ax^J +DW(y u )W(y ti )^ ti+1 , 

where x 2 .' a t . +i is the anti-symmetric part of the 2-tensor x 2 t . The first term only depends on 
Xj. t , and the second term can be simplified to 

2d 

J2 \W p ,W«}(y ti )xS£. 

p,q=l 

From the definition of (W 1 ) it is easy to see that the Lie bracket 

[W p , W q ] = Vp G {1, d} and Va 6 {d + 1, 2d} 

(and, therefore, it also vanishes for every p £ {d + l,...,2d} and q G {l,...,d} by antisymmetry). 
Each summand in (|3.10[) thus only depends on x 1 , x 2 and x 2 , , but not on the terms of x 2 corre- 
sponding to integrals between x} and x^; the same is hence true of the limit, y 8 j. 
We recall that y 2 1 is the limit as \D [s, t}\ — >• of 

[frA-n ® Vu,U + i + [W (y u ) 81 W (2/ ti )] x 2 )t(+l ] , (3.11) 

i:ti6D[s,t] 

so that in general y 2 does depend on the extension. However, by taking projections the dependence 
disappears. To see this just let fi : R e x R e -> R e denote the projection Pi (2/1,1/2) = Vi, so that 
y? = (Pi ® Pi) (y 2 ) ■ We then observe that 

(Pi ® Pi) (W (y t ) O (yt)) = (V (y t ) , 0) ® (V (y t ) , 0) , (3.12) 

and also the corresponding relation for y 2 . The claim then follows at once from p. lip and (|3.12l) . □ 

Corollary 3.2. Let v = X5 Xl + (1 — A) <5 X2 G M (P p .d)- There exists a unique fixed point /i G 
■Ml (Pp,e) of the map ^ v which is given explicitly by 

V = u r f* [X6 yi + (1 - A) 6 y2 ] , (3.13) 

where yi G WGQ, p (R 2<i ) , i = 1, 2 are the projections of the solution to i3. 6]) driven by any extension 
x of xi and X2 . 
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Proof. The previous proposition ensures (|3 . 1 3[) is well-defined. In other words, for every fixed 
realisation (y± (0) , y 2 (0)) of uf 2 the rough paths yi and y2 will not depend on the choice of 
extension. We muct check that this is the only fixed point. To do so, first note that the assumption 
on v implies that any fixed point must have the form 

where Zj, i = 1,2 are elements of GVL p (R e ) . Then, by the definition of the map ^> Vl Zi, i = 1,2 
must solve the RDEs 

dz, (t) = V° ( W (i)) drff + V ( Vi (t)) dx, (t) . 

Let x be any path in WG£l p (lR 2d ) whose projections are consistent with X! and x 2 . Then, since 
Zi and z 2 may both be written as solutions to RDEs driven by x, we may define in a canonical way 
(see [T3]) a path z in WGil p (M 2d ) , which has zi and z 2 as its projections, z is then the solution 
of the RDE driven along x . □ 

3.2. N-particle systems. We will later want to consider the propagation of choas phenomenon for 
rough differential equations, and this requires us to present the treatment of the previous subsection 
for a population of particles of arbitrary finite size N. We therefore suppose that the preference 
measure is now given by 

N 

v = \j6 Xi . 

i=l 

Analogously to the two-particle case (recall (|3.7[l ) we define vector fields W° and W = (W 1 , W Nd ) 
on R We = K e x .... x M e (this time as differential operators for notational ease) by writing y = 
(y 1 , y N ) e M Ne = R e x .... x Wt e and setting 

N e N „ 

^°(y) = EEE x ^ (y m > y l ) ^ ( 3 - 14 ) 

and, using the convention kd = d(mod<i) for k G Z, 

W> (y) = ±V^ (yi+LW-DZ-j) 1+[ g_ 1)/tiJ , i = l,...,Nd. (3.15) 
fc=i ^fc 

As before, we will be interested in rough paths in WGSl p (M wd ) whose projections contain each of 
the rough paths Xi, x 2 , xn, which together form the support of v. The following notation indexes 
the components of the extension in terms of the components of x 1; x 2 , ...,xjv. 

Notation 3. For each k inN and m = 0, ....,N — 1 define Ik,m; a subset of {I, ....,Nd} , by 

(ii, ik) G h, m iff {h, —, ik} C {md + 1, (m + 1) d} . 

We will let Ik denote the subset U^Zg/fe jm . 

We now formalise the precise sense in which {xi,x 2 , ...,xat} is related to any extension. 

Definition 3.3. If {xi,x 2 , ...,xjv} is a collection of rough paths in GQ P (W 1 ) , then we say that 
x in WGflp (R Nd ^j is a lift which is consistent with x 1? x 2 , x^y, if for every k = 1, |_pj its 
projections satisfy 

7r (<i,...,<*) (x)=x fci« 1 (mod«i),...,< fc (mod« l )) ) VC*!,....**) £ I k ,m, Vm = 1, N. 
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We now chose any lift x which is consistent with xi,X2, ...,Xjv ■ We want to show that if wc 
solve the RDE 

dy t = W° (y t ) dt + W (y t ) dx t> y (0) = ( Vl (0) , y N (0)) G R Ne , 

along x, then the output y will have the same projections irrespective of the initial choice of lift. To 
do so, we have to identify normal subgroup K of G n (R Nd \ so that {xi, X2, ...,xjy} can be identified 
with a path in the quotient group G n (R Nd ) /K. 

Lemma 3.4. For n G N let g" (R Nd ) denote the Lie algebra of G n (R Nd ) . Suppose t n (R Nd ) is 
the subset of Q n (R Nd ) defined by 

P (R Nd ) = [a e g" (R Nd ) : (e^a' 7 ') = 0,VJ e Ujf =1 / fc } , 

where, if I = (it, ik) , we write \I\ = k and e} := e* x (E> .... <8> e* fc . Let 

K n ^ R Nd^ . = exp (r (R^)). 

Then K n (R Nd ) is a connected Lie subgroup of G n (R Nd ) , I™ (R Nd ) is an ideal in g™ (R Nd ) and 
hence K n (R Nd ) is a normal subgroup of G n (R Nd ) . 

Proof. It is immediate that K n (R^) is a connected Lie subgroup. To prove that t n (R Nd ) is an 
ideal g" (R Nd ) we need to show that for any a in «" (R Nd ) and b in g" (R Nd ) wc have 

(e], [a,b} {1 ^ = for all I el k ,k = l,...,n. 
But this follows by noticing that 

(ej, (a ® 6) fe ) = J2 (eh a 1 ® ft*"') = £ (ej (0 , a') (ej (fc _ , 6*"') = 0. 
i=i z=i 

Where, for every / = 1, k — l,we have written / =: (/ (I) ,1 (k — £)) and used the fact that I G 

to deduce 1(1) G I; and L (k — I) G Ik-i- It is easily seen from this that (^j, [a, = 0. The 

assertion that if™ (R^J is normal then follows from the well-known correspondence between ideals 
of Lie algebras and normal subgroups of the Lie group (see, e.g., [30]). □ 

Remark 3.5. In a straight forward way we may uniquely identify any given collection of rough 
paths {xi,X2, ...,xjv} in GS1 P (R d ) with a path, which we denote by (xi, X2, xjv), in the quotient 
group: 

G lpi (RAM) /K lpi (R^rf) . 

(xi, X2, xjv) will then have finite p— variation with respect to the homogenous quotient norm (see 
|33j ). Any extension x £WG£l p (R Nd ) which is consistent with Xi, X2, xjv as described above, 
will then extend (xi, X2, xjy) in the obvious sense that 

^GLpJ (K Nd ),GLpJ (K Nd )/i?LpJ (E«d) (x) = (Xi, X 2 , Xjv) . 

We now prove the generalisation of Proposition 13.11 to the iV-particle system. 

Theorem 3.6. Let p > 1, y± (0) , ...,2/jv (0) G R e , and suppose that {xi,X2, ...,xtv} is a collection 
of rough paths in GQ p (R d ) . Assume that W° , defined by pTTJ\ ), and W = (W 1 , W Nd ) , defined 
by f3.15\) are, respectively, vector fields in Lip@ (R Ne ) and Lip 1 (R Ne ) for some /3 > 1 and 7 > p . 
For any q in [p, 7) let x be an element of WG£l q (R 2d ) which extends Xi, X2, Xjv in the sense of 
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iS.5)) . and suppose y be the unique solution in WGVl q {R Ne ^ to the RDE \3. 6)) driven by x. Then y 
has the property that its projections yi, ...,yjv to elements ofWGQ q (R e ) depend on Xi,X2, ...,xat, 
but not on the extension x.. 

Proof. From the LV Extension Theorem, there always exists an extension x of (xi,X2, ...,xjv) in 
WGQ q (R 2d ) for any q > p WGVl q (R Nd ) (and any q > p, if p is not an integer). Let us define an 
algebra homomorphism from the (truncated) tensor algebra tLpJ (R 7 ^) into the space of continuous 
differential operators, by taking the linear extension of 

F w (e^.ij = W 11 0...0W 1 ". 

Restricting F w to gLpJ (R^d) 

gives a Lie algebra homomorphism into the space of vector fields on 
M. Nd . An easy calculation confirms that 

te*(F w \ slpi{RNd) ) DtW (R Nd ) , (3.16) 

whereupon Theorem 20 of [33] shows that y 1 is independent of the extension of (xi,X2, ...,Xjv) to 
x. In general, y 2 , yL p J will still depend on the choice of lift. Nevertheless, the projections of y to 
the N paths yi,y2, ■•^yjv will be not do so. This is most easily seen when p S (2,3) by the same 
calculation as in (|3 . 1 1|) . □ 

Remark 3.7. Because each y, ; = &v°,v (f 1 . Hi (0) , x.j) , y,; solves an RDE driven by x^ G Gfl p (R d ) 
and the Universal Limit Theorem guarantees that y, ; is in fact an element of Gfl p (R e ) . This ob- 
servation will be useful later on. It follows from Theorem Iff. 61 together with a suitable elaboration 
of the arguments of Corollary \3. 2\ that 

N 

1=1 

is the unique fixed point of <lv 

4. A FIXED-POINT AND CONTINUITY THEOREM 

We now want to consider the case where the preference measure v is a non-discrete measure on 
rough path space. The main problem we address is to find a condition on v to force the existence 
of a unique fixed point to the map \£ r 1/ . A key feature will be the use of estimates controlling: 

Pp~var;[0,T] (y\ y 2 ) i 

the /9 p _„ ar; [ 0! T]-distance between two RDE solutions y 1 and y 2 driven by x. These estimates need 
have two properties: they need to be Lipschitz in the defining data (starting point, vector fields 
etc) and the Lipschitz constant must have integrable dependence on x, when x is realised according 
to a wide class of measures. Classical RDE estimates satisfying the first of these criteria, the latter 
needs more work. For example in [20] the authors have proved estimates of the form 

Pp-var-%T] (y\y 2 ) < (*) ex p(ll X llp-™r;[0/r]) , C 4 - 1 ) 

where the terms (*) incorporate the data. The drawback of this estimates is that the right hand side 
fails to be integrable, for example when x is the lift of a wide class of common process including 
Brownian motion and fractional Brownian motion with H < 1/2. Fortunately, it it possible to 
replace |x| |p_^ ar .[ T ] in (|4.1[) by a quantity called the accumulated a— local p-variation (sec below). 
By then making use of the recent tail estimates in [5], we are able to cover these interesting examples. 
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4.1. Lipschitz-continuity for RDEs with drift. We recall the definition of the following func- 
tion from [5]: 

Definition 4.1. Let a > and I C R be a compact interval. Suppose that ui : I x I — > R + is a 
control. We define the accumulated a—local ui-variation by 

M a ,;(w)= sup ^2 ^ (U,ti+i) ■ 



The following lemma is a Lipschitz estimate on the RDE solution (with drift), when we vary the 
defining data of the differential equation. 



Lemma 4.2. Let 7 > p > 1 and (3 > 1. Suppose x is a weakly geometric p— rough path in 
WGflp (M- d ) , and assume that 7 1 and j 2 are two paths which take values in some Banach space E, 

40, T] 



and belong to C 1 ""^ ([0, T], E) . Then u : A [0T] -> M+ de/med by 



2 

l — var:[s } ■/ I I 



( s ; *) : ~ X! 1 1 7 II l-«or;[*,tl + " 



is a control. Furthermore, if V = (V 1 , is a collection of vector fields in Lip 1 (R e ) , and V° 
is in LipP (K e , L (E, M e )) , then for i =1,2 tfte i?D£s 

dyj = V {yi) d*l + V° (yf) drf t) 

my* = wo 

/lave unique solutions. And for every a > and some C = C(w,a) > 0, we also have the following 
Lipschitz-continuity of the solutions: 

Pp,u (y\y 2 ) < C [\yl -yl\ +p hu (l\l 2 )+P P ,u (x\x 2 )] exp (CM Q , [0 , T] (w)) . 

Proof. The proof is obtained by following the arguments of Theorem 12.10 of [30J on RDEs with 
drift; two enhancements are necessary. The first is allow the drift term to take values in an arbitrary 
(infinite dimensional) Banach space. This is elementary, because in the current lemma 7 1 and 7 2 
have bounded variation, and hence classical ode estimates can be used everywhere. The second, 
more subtle, enhancement is to end up with the accumulated a— local w-variation featuring in the 
exponential on the right hand side (as opposed to the usual uj (0,T)). For this we refer to [3] and 
Remark 10.64 of [30]. □ 

By exploiting the relationship between p p ^ and p p - uar we can obtain a Lipschitz estimate in 
Pp-uar-distance: 



Corollary 4.3. With the notation of, and under the same conditions as, Lemma \4.2\ we have 

Pp-var,[0,T\ (y\ y 2 ) 

< Cu(0,Tf [\yl-yl\ + | |V - 7 2 || 1 _ wr;[0)T] + | (x 1 - x 2 | | p _„ ar;[M] ] exp (CM a , [0 , T] (w)) , 



for some N > 0. 
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4.2. Measure-valued paths. For the current application, the main interest in these Lipschitz 
estimates will occur when the space of probability measures M. (S Pje ) is embedded in a Banach space 
E. In the typically case 7 will then be constructed from \i G M (P p . e ) by setting 74 := f Q p s ds. For 
the moment, we develop this more abstractly by letting Lip 1 (S*)*denote the dual of Lip-1 functions 
(that is, the bounded Lipschitz functions) on a metric space (S, d) . There is a canonical injection 
p t-» from Ai (S) into Lip 1 (S)* defined by the integration of functions in BL (S) against p : 

r M (^) = (T MJ ^) := / <f>(s)p(ds). (4.2) 
Js 

In this setting, two metric spaces will be of special interest as already mentioned in Section [2] 
The first is the step-iV free nilpotent group with e generators, G N (M e ), with the (inhomogeneous) 
metric it inherits from the tensor algebra: 

rfjv(g,h):= max |7r,(g-h)|. 

i— 1,...,N 

The second is the space of geometric p— rough paths Gf2 p (R d ) equipped with p p - va r:[o.T] ■ 

The next lemma examines the regularity of the paths which result from (|4.3[) the pushforward 
of p under the evaluation maps, i.e. 



Vt = C0t)* A* e M (S p , e ) . 



(4.3) 



Lemma 4.4. Suppose p > 1 and iei fee an element of M.\ (P p , e ) ■ Let (Mt)tg[or] ^ e ^ e 

m ./VI {Sp. e ) defined by h4-3\ ), and C^ t )te[oT] ^ e pas^i Ml Lip 1 {S p . e ) obtained by the injection of 

A4 (S Pl e) into Lip 1 (S Pye )* . Then for every < s < t < T we have 



(Ys, Yt) P (dy) < / Pp-var;[s,t] ih y) M (<t>0 j 



where y u — ip u (y) . in particular, t 1— > T Mf is a continuous path in Lip 1 (S PtC ) . 

Proof. Take £ Lip 1 (S p ) with |0| | £ip i( S \ = 1. The result then follows from the proceeding 
calculation: 



1(7^-7^,0)1 



0(g) [fa)./*] {ds) 



0(g) [(V-).A*] (*) 



< 



< 



° ipt - & ° ips) (y) v (dy) 
d[ P \ {i> s (y),*l>t (y))/i(dy) 
Pp— Bor;[s, *] (i,y)M(dy) 



The right hand side is finite since p is in A^i (P p , e ) , and (by the dominated convergence theorem) 
it tends to zero as \t — s\ tends to zero. □ 



It follows from this lemma that 

7s, t 



■= / T ti r dr ■= / l[ S) t] (r)T Mr dr G Lip 1 (S p , e )* 



(4.4) 
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is well-defined for every (s,t) C [0,T] , where the integral is understood in the sense of Bochner 
integration. For any <f> S Lip 1 (S Pje ) , standard properties of the integral yield that 



( 7s , t ,0>=/ (T fir ,<t>)dr= / 0(g)fi r (dg)dr. (4.5) 

The space P p%e carries with it an implicit time interval [0,T] which we suppress in the notation. 
Occasionally, we might want to make this explicit by writing P Pi e,T- For example, if we start with a 
probability measure [i in M. (P^.t) we will need to consider its restriction, fi\[o.t] to a probability 
measure in A4 (P p , e ,t)- We then let Wt denote the Wasserstein metric on M. (P p , e ,t), and write 
W t (m\m 2 ) to mean W t (n\o,t], M 2 |[0,t])- 

Corollary 4.5. Suppose p > 1 and let T > 0. Assume pi 1 and ji 2 are two elements of M.\ (P p .c,t) , 
and for i = 1, 2 let j l : A[ .t] — > Lip 1 {S p , e )* be the function defined by 

iU = l T ^ dr - 

Then for every < s < t < T we have that 

H,t-llt\ Lip , (Sp ^ <C{t-s)W t ^\f). (4.6) 

In particular, 

IK-^IL^o.t] <cJ o T W t (v\^)dt. (4.7) 
Proof. Let <fi G Lip 1 (S P:t ,) with H^H^^^ y = 1, then from (|4.5p we can deduce that 

|<7i,t-7«,t>^)| = 



0(g) M J(dg)- / ^(g) Mr (dg) 
<(£-s)sup / d LpJ (g\g 2 ) 7r r (dgSdg 

re[s,t] JSp. e XSp e 



where 7r r is an element of .Mi (5p, e x S p ^ e ) with marginal distributions pi 1 and /i 2 . For any such 
7r r and every re [s, £] we have 

/ d [pi (g 1 ,g 2 )n r (dg 1 ,dg 2 ) <C P P -var:[o,t](y 1 ,y 2 )^(dy 1 : dy 2 ) , 

JS p , e XS p ,e J Pp, c ,tXP P ,e.t 

where ir in A'fi (Pp, e ,t X Pp,e,t) is any coupling of /Lt | ro,t] an( i M 2 |[o,*]- This yields 

SUP ( d M (g 1 , g 2 ) 7T r (dg 1 ,dg 2 )< W t (fl 1 , M 2 ) , 

r(E[s,t] JS p , c xSp, e 

which implies (|4. 6[) at once. Deducing (|4.7[) from (|4.6[) is then elementary. □ 

4.3. A fixed-point theorem. Suppose that p > 1 and x is an element of G£l p (R d ) , then we write 
w x for the control induced by x via 

MM) = \\x\\ P p - va r;[s,t] ■ 

The following lemma gives a useful way of controlling cj x (0, T) in terms of the a-local p-variation. 
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Lemma 4.6. For any x in Gil p (K rf ) and any a > 0, we have that 

M P p -var;lo,T] = ^ x (0, T) < 2"~ V max { 1 , a" p M Q , [0 , T ] K) P } 

Proof. Fix a > 0, and let D — (ti : i — 0, l....,n) be an arbitrary partition of [0,T] . We aim to 
estimate 

n n 

^2 || x *i-i,t,|| P : = ^2 d cc (xt j _ 1 ,x ti ) p 
1=1 1=1 
Let and be any two consecutive points in D, and define <to — tj—i ancl 

cr i+1 = inf {t > ai : uj (ai, t) = a} A £j 

for i G N. Define 

Nccfa-utj] ( w x) = sup {n G N U {0} : a n < tj} 
A simple calculation shows that fjv = tj if N = iV^, .^.j (w x ) + 1, and therefore 
/w+i \" 

iK-i^ir - ( s ii xct «-i^ii) 

JV+l iV+1 

< (jv + i)^ 1 ^ i |x ffi _ l!CTj i i p < (7v a>[ o, T ] k) + i) p_i i ki-i.^ i r • 
i=i i=i 

Using this observation it is easy to deduce that 

n 

I | X *i-i>*« I T ^ (^.[O.T] (Wx) + 1)^ M a ,[o,T] K) • (4.8) 

1=1 

The claimed bounded follows by first noticing that N a jo,t] ( w x) < Qf _ Af« [o,t] i ail d then taking 
the supremum over all partitions D in (|4.8p . □ 

In order to prove the fixed point theorem we require integrability on the preference measure. 
The subset of M (P p ,d) for which the fixed-point theorem will hold is described by the following 
condition. 

Condition 1. Let p > 1. v will denote a probability measure in A4 (P p .d) > ® n d 4>v will be pushfor- 
ward measure in Ai ([0, oo)) defined by 

<f> v := [M h[0 ,T\ (w-)]. ( v ) ■ 
We will assume that <f> v has well-defined moment-generating function; i.e. for every 9 in R we have 

exp [By] <j> v (dy) = / exp [flMijo.r] (w x )] v (dx) < oo. 

[o,oo) Jp Pid 

Remark 4.7. If Condition[J\is in force, then := [M q [o,t] ty) a ^ so have a well-defined 

moment generating function for any a in (0, 1). 

For the reader's convenience, we recall some notation from Scction[3] "J : M.\ (P p , e ) — > -Mi (Pp,e) 
is defined by 

* = : /i ^ [@v°,v (A*, 0] * («o xv)eMi (Pp, e ) , 
and fixed points of ^ v correspond to solutions of the nonlinear McKcan-Vlasov RDE 

dY t = V (Yn dX t + V° (Yf) drf (4 9) 

Law(Y) = /i, Law(Yo) = u a 



EVOLVING COMMUNITIES 



17 



We now formulate and prove our main existence and uniqueness theorem for solutions to (|4.9[) . 

Theorem 4.8. Let 7 > p > 1 and f3 > 1. Suppose v be an element of Mi (P p ,d) which satisfies 
Condition^ Let {Yq,~K.) be a random variable on a probability space (tt, J 7 , P) , taking values in 
R e x P Pt d and having law Uq x v. Then for any collection of vector fields V = (V 1 ,..., V ) in 
Lip~< (K e ) and V° in Lip? (R e , L (Lip 1 (S p . e )* , M e )) , there exists a unique solution to the nonlinear 
McKean- Vlasov RDE fTl). 



Proof. The space (P p , e , Pp) is complete, and hence (see, e.g., [H]) so is (Mi (P p , e ) ,W) . Suppose 
fii and H2 are in Mi (P P , e ) , and let 

y 1 = e V o iV ( M i,y ,x). 

Using Corollary 14.31 together with Lemma |4~51 we obtain for any a in (0, 1] the bound 

P P -var; [Q ,T] (Y 1 , Y 2 ) < C | \^ - 7 W I ll-„ar ; [0,T] CXp (CM.^, (w X )) , 

for some C = Ci (a) > 0. It follows from Corollary 14.51 that 

ll7 Ml -7 M2 lli-„ or;[ o,T]^ [ W t (m,fi 2 )dt. 

•J 

And therefore by taking expectations in the previous inequality we obtain 

WT(* v (»l),VA»2))< E [p P -var;[0,T] (Y\Y 2 )] 

i-T 

w t (m,n2) dt, 



< CE [exp (CM a ,| 0iT] K))] f 

Jo 



where Condition [T] ensures that the right hand side is finite. It from a standard Banach-type 
contraction argument that the map $ v has a unique fixed point. □ 

In light of the conditions of this theorem, it is useful to make some observation about the type 
of processes which satisfy the key integrability condition (Condition [T]). In the recent paper [5] we 
consider a continuous Gaussian process X = (X 1 , ..,X d ^ with i.i.d. components such that: 

(1) X has a natural lift to a geometric p-rough path X; 

(2) The Cameron-Martin space associated to X has the embedding property 

U C«- var ([0,T] ,lR d ) 

for some l/p+l/q> 1. 
We then prove that for some ry > we have 

exp 7yM QjM (w x ) 2/9 



E 



< 00. 



This class of examples is rich enough to include fractional Brownian motion H > 1/4 (for which q 
can be chosen to ensure 2/q > 1), and other examples of Gaussian processes which are genuinely 
rougher than Brownian motion (see |19j). The importance of the Lipschitz estimate in Corollary 
14.31 can now be grasped more clearly. Since, as an immediate corollary, we see that Condition [T] 
holds for the class of measures described. 

Remark 4.9. In some recent work [2], a flow-based approach is used to derive continuity estimates 
for RDEs. An existence and uniqueness theorem is proved, under the following condition on X; for 
some family of random variables {C s :s£[0,T]}, which is bounded in L 1 , 

\E[X k a>t \r a ]\<C B (t-8),V[8,t]C[0,T\, fc = l,...,LpJ. (4.10) 
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This requirement forces some structure upon the sample paths of X , (for example: smoothness, or 
independence of increments). It does not hold in general for the examples illustrated above, where 
the sample paths are less regular than Brownian motion. Indeed if X is fBm with H < 1/2, then 
we have 

r x E [Xlt] = t 2 "- 1 t oo as 1 1 0, 

which violates j-j-lO^ when s = 0. By contrast, the exponential integrability required in Condition]]^ 
holds both in this example, and for the much wider class of Gaussian processes highlighted above. 

4.4. Continuity in v. Suppose we have a set of preference measures and for each measure in the 
set the conditions of Theorem 14.81 hold, so that ^> v (•) has a unique fixed point. A very natural 
question is to ask about the stability properties of this map. The rough path setup is well-suited 
to tackle this sort of problem. To this end, let K : (0, oo) — > (0, oo) be a monotone increasing 
real-valued function. Define a subset of Mi {P p ,d) by 

£ (K) = {veMi (P M ) : Vfl e (0, oo) , </v (6) < K (8)} , 

where, as above, 

VV (0) = f cxp (By) cj> v (dy) . 

7[0,oo) 

It is easy to see that £ (K) is a closed subset of Mi (P P ,d) m the topology of weak convergence of 
measures. 

Lemma 4.10. Let K : (0, oo) — > (0,oo) be a monotone increasing, then the map 

3 : £ (K) M\ (P e ,d) 

3 '. v i — \ fixed point of ^ v (•) 

is well-defined and continuous in the topology of weak convergence of measures on £ (K) . 

Proof. Suppose that is a sequence of measures in £ (K) such that v n v € £ (K) as 

n — > oo; we will show that Wt (3 (v n ) , 3 (v)) — >• as n — > oo. By Skorohod's lemma there exists a 
probability space carrying: (i) an M e -valucd random variable Yq with law uq, and (ii) a sequence 
of (P Pi( z- valued) random variables (X 1/ ")^° =1 and X", such that X" has law v, X" n has law v n for 
every n, and 

X"" — > X" a.s. in p p - V£LT - 

Let 3(i/„) = fi n , 3(v)=n and Y"» = e V o jV {fi n , Y , X v »), Y" = Q v oy (ft, Y , X") . We then 
have that 

Law(Y I/ ") = /i„ = [6y ,v (/•*«> •>•)]* ( u o,^n) , 
and similarly for Law(Y") . The estimates of Corollary 14.31 and Lemma 14.61 show that for some 
non-random C\ > : 

p(Y Un ,Y V ) < & [\\>f« ^Wl-varW + ll X '" - XU W P -var ;M ] «P (ClM a>[0 , T] , 

where 

w^r*,*) = iix^ip r ,, + iixt r .,. 

V ' / II iip—var;[s,t\ II Ijj— war; [s,t\ 

Taking expectations and then making use of Corollary 14 . 5 1 and Condition [1] it is easy to derive that 

W T (fi n ,fi) <C 2 [ W t (fi n , fi) dt + da n , (4.11) 
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where 



a„ := E 



IX 1 "* - X"|| p _ rar;[0iT] exp (CM aA0>T] (oA'"))_ . 

A simple argument using the definitioifl of £ (K) shows that a n — > as n — > oo. Finally, we use 
Gronwall's inequality in (|4.11[) to give 

Wj 1 (fx n , n) < C^a n exp (C3T) — > as n — > 00. 



□ 



5. Applications 



As an application of our uniqueness theorem, we prove the classical propagation of chaos phe- 
nomenon (see Sznitman |42j ) for the finite interacting particle system. This is the observation that, 
granted sufficient symmetry to the interaction and initial configuration, then in the large-population 
limit any finite sub-collection of particles resembles the evolution independent particles, each hav- 
ing the law of the nonlinear McKean-Vlasov RDE. To make progress, let v be a fixed preference 
measure in A^i (P P ,d) which satisfies Condition!]] Assume that {(X%Yq) : * £ N} and (X,Yb) arc 
i.i.d. R e x P p ,d— valued random variables each with law uq x v and defined on the same probability 
space (O, J 7 , P). Suppose pi is the unique fixed point of ^ v (•) , and let 

Y = 9^ y (n,Y ,X). 

From Lemma 13.21 we can interpret the trajectories of individual particles in the community 

1 N 

dY * N = Jj E CT Y ?' N ) dt + V W) Y o N = Y o 

3=1 

as projections of the solution to a system of rough differential equation driven by any rough path 
in P Pt Nd which consistently lifts X 1 , ...,X . Indeed, we showed that Y l is then well-defined as 

Y^ N = @v°,v (/Ar \x% i = l,....,N 

where V° in Lip 13 (R e , L (Lip 1 (S Pye )* , R e )) is defined in terms of the interaction kernel a by 

V° (y) = a (y,-)), (5.1) 

and /j, N — [i N (uj) is the empirical measure 

JV 



For every N, {Y 1,N , Y 2,JV , Y'*'"} is an exchangeable system of random variables, and a classical 
result of [42| shows that propagation of chaos is equivalent to proving that 

Remark 5.1. Explicitly, this assertion says that the law of the random variable 

fi N (u)eMi(P P ), 

which is a probability measure in Ai\ (Mi (P p )) , converges weakly as N — > 00 to the probability 
measure 8^, which is the law of the constant random variable fi. 



In particular boundedness (uniform in n) of the moment generating functions of u n 
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Theorem 5.2. Let j > p > 1, (3 > 1 and j/o S R e - Suppose that is is a given preference measure 
in M.\ (Pp.d) which satisfies Condition^ Assume that the vector fields V = (V 1 , V d ^j belong to 
Lip* 1 (R e ) , and V° defined by i f 5. 1\) is in Lip? (R e , L (Lip (S p , e )* ,W)) . Let pi denote the unique 
fixed point of map $ '„ (•) which results from Theorem \4-8\ Assume further that {(Yq, X 1 ) : i s N} 
is a collection of i.i.d. R e x P p ,d— valued random variables, with law uo x v 1 defined on the same 
probability space. For each N £ N let |Y l ' W : 1 = 1, ....,iV} be the solution to the particle system 

1 N 

^ N = ^J2 (T Y t' N ) dt + v ( Y t) Y' d N = y„«el e . 
i=i 

Then as N — >■ oo 

1 N 

=•■ f N => * ( 5 - 2 ) 

j'=i 

and the particle system exhibits propagation of chaos. 



Proof. Fix N g N and suppose i e {1, ...N}. Let Y* = Qv°,v (m>^ojX*) so that the law Y* is /i, 
and {Y 1 , Y w } are AT independent copies of the solution of the rough McKean-Vlasov equation. 
Then using Corollary 14. 31 and Lemma T4. 61 we have (for a in (0, 1]) that 



Pp-var;[0,T] (Y*, Y l,N ) < C\ 

Using Corollary 14.51 we observe 



l-var;[0,T] 



exp (CiM Qi[0 ,r] (wxO) • ( 5 - 3 ) 



<C 2 I W t (fi,n N )dt 



l-var;[0,T] 

and hence by summing (|5.3p over i = 1, JV we obtain 

V/l rnm, (Y*. Y i,JV> l < ^ / 

TV 

Let 



o 



^Epp-™ r ;[o,T] (Y*,Y*-*) < C 2 / Wi ( / ,, M A, )^l^exp(C 1 M Qi[0jT] (a*,)) 

i=l *' i=l 



N 



Then we also have 



W T U, fj, N ) < W T U H n ) + W T (fl N , fi N ) 



And, on the other hand using (|5.4p we have the bound 

T N 

W T (fl N , M W ) < C 2 [ W t (/x, M W ) dt± J2 ex P ( C ^[o,T] M) 

*' i=l 

Putting (|5.6j) into (|5.5|) and using Gronwall's lemma we deduce that 



W T (h, I* n ) < C 2 W T (h, fl N ) exp 



C 

Tf H CX P ( C 'l M o,[0,T] OxO) 



N 



(5.4) 



(5.5) 
(5.6) 
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It is a simple matter to conclude from the strong law of large numbers that both Wt (m, — > 
a.s., and 

1 N 

— J^exp (CiM Qj [ ,t] OxO) -> # [exp (CiM Q , [0jT] (u x ))] < oo, 

i=l 

a.s. as iV — > oo. It is then easy to deduce (|5.2|) . Propagation of chaos is then a consequence of the 
classical result of [55] we cited earlier. □ 

There are a number of follow-up results that seem worth pursuing. For example, Sanov-type 
theorems a la Dawson-Gartner [3] will be possible for (|5.2p in the weakly interacting case. Indeed, 
the presence of the rough path topology, in which the universal limit theorem guarantees the 
continuity of the Ito map seems to simplify things greatly. We will return to these discussion in 
future work. 
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